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Abstract. Recent models of rotating massive stars including magnetic fields prove it difficult for the cores of single stars to 
retain enough angular momentum to produce a collapsar and gamma-ray burst. At low metallicity, even very massive stars may 
retain a massive hydrogen envelope due to the weakness of the stellar winds, posing an additional obstacle to the collapsar 
model. Here, we consider the evolution of massive, magnetic stars where rapid rotation induces almost chemically homoge- 
neous evolution. We find that in this case, the requirements of the collapsar model are rather easily fulfilled if the metallicity 
sufficiently small: 1) Rapidly rotating helium stars are formed without the need to remove the hydrogen envelope, avoiding 
mass-loss induced spin-down. 2) Angular momentum transport from the helium core to hydrogen envelope by magnetic torques 
is insignificant. We demonstrate this by calculating evolutionary models of massive stars with various metallicities, and derive 
an upper metallicity limit for this scenario based on currently proposed mass loss rates. Our models also suggest the existence 
of a lower CO-core mass limit of about 10 M Q — which relates to an initial mass of only about 20 M Q within our scenario — for 
GRB production. We argue that the relative importance of the considered GRB progenitor channel, compared to any channel 
related to binary stars, may increase with decreasing metallicity, and that this channel might be the major path to GRBs from 
first stars. 
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1. Introduction 

As long gamma-ray bursts (GRBs) are observable, in princi- 
ple, from redshifts as high as 20 (e.g. Wijers et al.Q]598 1, their 
study may yield break-throughs in our understanding of the 
early universe. Their association with the deaths of massive 
stars (e.g. Hjorth et al. 2003 Stanek et al. 2003) renders them 
as tracer of the cosmic history of star formation (see Bromm 
& Larson 2004 for a review). Their possible use as standard 
candles (e.g. Friedman & Bloom 2005 , Lamb et al. 2005 ) may 
turn the into the farthest reaching cosmic yardstick. However, 
our present lack of understanding the progenitor evolution of 
long GRBs hampers their deeper understanding and quantita- 
tive exploitation for cosmology. 

According to the currently favoured collapsar scenario, rel- 
ativistic jets are formed when rapidly rotating iron cores of 
massive helium stars collapse to form black holes and accre- 
tion disks around them (Woosley 1993 ). However, recent stellar 
evolution models which include angular momentum transport 
from the core to the hydrogen envelope by magnetic torques 
(Spruit 2002) indicate the most single stars end up with too 
slowly rotating cores. Especially during the giant stage, where 



strong differential rotation persists between helium core and 
hydrogen envelope, angular momentum is efficiently removed 
from the core (Maeder & Meynet 120031 Heger et al. 120051 
Petrovic et al. 120051 see also below). Although some non- 
magnetic stellar models are found to fulfill the necessary con- 
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ditions for collapsar (Heger et al. 2000; Petrovic et al. 2005, 
Hirschi et al. 2005), they predict a much too high GRB vs. 
Type Ibc supernova ratio (e.g., van Putten 2004). The inclu- 
sion of magnetic torques is further required to reproduce the 
observed low rotation rates of young neutron stars (Heger et 
al. 2005) and white dwarfs (Suijs et al. 2005). 

The formation of collapsars is further complicated by con- 
sidering mass loss. For rather large mass loss rates, a rotat- 
ing star is efficiently spun-down even without magnetic effects 
(Langer 1998 ). This is in particular true also for massive helium 
stars, which, in our Galaxy, appear as Wolf-Rayet star with 
wind dominated spectra: any long-lasting Wolf-Rayet phase at 
solar metallicity will likely lead to smaller spins than required 
by the collapsar model. However, Wolf-Rayet mass loss rates 
may be smaller for smaller metallicity (Nugis & Lamers 2000 
(NL00); Crowther et al. 120021 Vink & de Koter 120051 . The 
problem then is that the winds of potential Wolf-Rayet star 
progenitors also weaken for lower metallicity (Kudritzki et 
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al. 119891 Kudritzki 2002 1, with the consequence that the Wolf- 
Rayet stage may actually never be reached (Heger et al. 2003 
Marigo et al. l2"U0H Meynet & Maeder l^UOH . 

Here, we explore an evolutionary channel which can avoid 
both, the problem of magnetic torques, and the mass loss 
problem, at low metallicity: chemically homogeneous evolu- 
tion. Maeder ( 1987 1 found that if rotationally induced chemi- 
cal mixing in massive main sequence stars occurs faster than 
the built-up of chemical gradients due to nuclear fusion, the 
chemical gradients may remain small throughout core hydro- 
gen burning, transforming the main sequence star smoothly 
into a helium star of the same mass. Even though the physics 
of rotationally induced mixing has been refined since then, 
quasi-chemically homogeneous evolution is still found if very 
rapid rotation is adopted (Langer 1992 Heger & Langer 2000 
Maeder & Meynet 2000 1, and remains possible also when mag- 
netic torques are included (see below). 

In quasi-chemically homogeneously evolving models, the 
helium star is born without, or with only a small hydrogen- 
rich envelope. Magnetic torques can therefore remove only lit- 
tle angular momentum from the helium star. And at low enough 
metallicity, the mechanical mass loss induced spin-down may 
be avoided due to the weakness of the helium star winds. After 
we explain our numerical methods and physical assumptions 
in Sect. |3 we present models of the homogeneous evolution of 
rapidly rotating massive stars at low metallicity with the effect 
of magnetic torques included in Sect. [3] We discuss the poten- 
tial implications of our results for an understanding of GRB 
progenitors in in Sect.|4] 

2. Methods and physical assumptions 

Stellar models are calculated with a hydrodynamic stellar evo- 
lution code (cf. Petrovic et al. 2005 , and references therein), 
which includes the effect of the centrifugal force on the stellar 
structure, chemical mixing and transport of angular momen- 
tum due to rotationally induced hydrodynamic instabilities, and 
the transport of angular momentum due to magnetic torque 
(Spruit 2002; Maeder & Mevnet l2005l . 

Stellar wind mass loss for O stars is calculated accord- 
ing to Kudritzki et al. ( 1989). Wind mass loss rates of helium 
rich stars (WR stars) are currently very uncertain. Recent WR 
star wind models, which include the effect of wind inhomo- 
geneities (clumping), give significantly lower values than ear- 
lier estimates, e.g. by Hamann et al. ( 1995 hereafter HWK95), 
as shown in Fig.^ Here, we use two different Wolf-Rayet loss 
rates (WR1 and WR2), as follows: 



WR1 : log^/fMoyr 1 ]) 



and 
WR2 



log^/tMoyr- 1 ]) 



-12.73 + 1.5 log L/L Q -2.85X S 
+O.51og(Z/Z ), (1) 

-13.13+ 1.5 logL/L -2.85X s 
+0.861og(Z/Z G ), (2) 



where X s denotes the surface hydrogen abundance. Both rates 
correspond to HWK95, but reduced by factors 6 and 15, and 
the metallicity dependence is M a Z 5 and Z 86 for WR1 
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Fig. 1. Wolf-Rayet wind mass loss rates as a function of the 
stellar luminosity for a surface composition of X = 0.15, 
Y = 0.83, and Z = 0.02. The mass loss rates by Hamann et al. 



( 1995 , HKW95) and Nugis & Lamers (2000, NLOO) are given 
by the solid and dashed lines, respectively. The upper and lower 
dotted lines denote the HKW95 rates divided by a factor of 6 
and 15 respectively, as assumed in this study. The filled circle 
and the open triangle are the mass loss rates given by Vink & 
de Koter ( 2005 1 for WN stars, with two different assumed val- 
ues of p, which determines the wind velocity law (see Vink & 
de Koter 2005 for the exact definition of /3). 

and WR2, respectively, respectively (see Fig. [Q. The former 



is commonly used in stellar models (e.g. Heger et al. 2003 
Hirschi et al. 2005 1, while the latter is based on the recent WN 
stars wind models by Vink & de Koter ( 2005 1. 

The enhancement of stellar mass loss when the star ap- 



proaches the Q-limit is considered following Langer ( 1997 



M = M(V l0t = 0) 
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1 - V/Vcrit 
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where 



V a 



GM 



-(i-D; r 



kL 



(3) 



(4) 



R v _ " " 4ncGM 

We include the effect of the enrichment of CNO elements at the 
stellar surface on the opacity, and its influence on the Eddington 
factor F. 

Initial rotation rates of our models are chosen such that ro- 
tational velocity at the equatorial surface has a specified frac- 
tion of the Keplerian value: /k = Vinit/Vkepier- We compute 
models forZ = 10~ 5 , Z = 10~ 3 , and one model with z - 0.02 
(see Table Q. The chemical composition is assumed to have 
the same proportionality with those of solar metallicity. More 
details about the numerical code can be found in Petrovic et 
al. (2005 1 and references therein. 



3. Results 

Table Ogives an overview of the computed model sequences 
and their basic properties. All model sequences with Z = 10~ 5 
and Z = 0.001 are followed until the exhaustion of central car- 
bon or the beginning of neon burning, after which the evolu- 
tionary times are too short (< 1 yr) to meaningfully change the 
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Table 1. Basic model properties: Sequence number, M^: initial mass, Z: metallicity, M WR : option for the WR wind mass loss 
rate, f K : initial rotational velocity in units of the Keplerian value, Vj n i t / V C rit : the ratio of the initial rotational velocity to the critical 
velocity at the equator (Eq. 4), F s ,H_ e nd : surface helium abundance at core hydrogen exhaustion, Msn^. final mass, Mqo- CO core 
mass at core carbon exhaustion, Ji nlt : initial angular momentum, /final : final angular momentum, 73m : mean specific angular 
momentum in the innermost 3 M Q , Xisj.max: maximum nitrogen surface mass fraction during the evolution. 
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Fig. 2. Evolutionary tracks for models of Seq. 3 (dashed line) 
and Seq. 4 (solid line), in the HR diagram (cf. Table 1), from 
the zero age main sequence to core carbon exhaustion. 



quantities given in Table ^ until collapse. The calculation for 
the model sequence B2 (Z = 0.02) is stopped at the end of the 
main sequence. 

Fig.EJshows that models of Seq. A3 (M = 40M o ,/ K = 0.2) 
evolve in a similar way to non-rotating models of 40 M and 
Z - 10~ 5 by Meynet & Maeder (2002 1, rather than their rotat- 
ing models which evolve further redwards. This is because in 
our models, the degree of differential rotation is much weaker 
due to magnetic torque. In Seq. A3, helium enrichment at the 
end of main sequence (AY = 0.03) is somewhat larger than in 
the 40 M rotating model with V init = 300 km s -1 of Meynet 
& Maeder (AY = 0.02). Recent solar metallicity models with 



magnetic fields by Maeder & Meynet (2005 1 show larger he- 
lium enrichment than their non-magnetic rotating models. Note 
also that our adopted initial helium abundance is slightly higher 
(Y = 0.24) than theirs (Y = 0.23). 
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Fig. 3. Evolution of the internal structure of Seq. A4 from zero- 
age main sequence to the neon burning. Convective layers are 
hatched. Semi-convective layers are marked by dots. The gray 
shading gives nuclear energy generation rates in log scale, as 
indicated on the right side. The topmost solid line denotes the 
surface of the star. 



On the other hand, models of Seq. A4 - which initially ro- 
tate 2 times faster than in Seq. A3 - evolve bluewards due to 
the strong helium enrichment of the envelope via rotationally 
induced chemical mixing. The evolution of the internal struc- 
ture of the star is shown in Fig. [3] The stellar mass decreases 
rapidly when the star enters the WR stage. At core hydrogen 
and core helium exhaustion, the stellar contraction leads the 
stellar surface to reach the Q limit, inducing strong centrifu- 
gally supported mass loss. In the end, the CO core is nearly 
four times more massive than in Seq. A3, even though the star 
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Fig. 4. upper panel: Mass fraction of chemical elements after 
core carbon exhaustion in Seq. A3, as a function of the mass 
coordinate, lower panel: Same as in the upper panel, but for 
model Seq. A4. 
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in Seq. A4 loses more mass. Fig. 0] compares the chemical 
profiles of Seq. A3 and Seq. A4 at the central carbon exhaus- 
tion. In Seq. A4, only a small amount of hydrogen is left at 
the stellar surface, but significant amounts of primary nitrogen 
(Xi4 N =i 10~ 2 ) are present. In Seq. A3, the nitrogen abundance 
at the surface is about 410~ 6 at this stage. 

We note that in Seq. B2 (Z = 0.02 and f K = 0.5), the he- 
lium enrichment is less than in the corresponding low metal- 
licity models (see Table [Q. The reason is that chemical mix- 
ing becomes inefficient as the star is significantly spun down 
via strong stellar winds, given that the mixing time scale is 
inversely proportional to the stellar mass and spin rate (e.g. 
Maeder & Meynet 2000 1. This results indicate that homoge- 
neous evolution is favored at lower metallicity, since the stars 
suffer less mass loss. 

The chemical mixing affects the internal transport of an- 
gular momentum in our models. Fig.|5]shows that much more 
angular momentum is removed from the core in Seq. A3 than in 
Seq. A4. In Seq. A3, while stronger chemical gradient between 



Fig. 5. Mean specific angular momentum over the shells as a 
function of the mass coordinate at 4 different epochs in Seq. A3 
(upper panel) and Seq. A4 (lower panel). Note that the equato- 
rial specific angular momentum is larger by a factor of 1 .5. The 
thin two-dotted-dashed line denotes the angular momentum 
of the last stable orbit around a non-rotating Schwardzschild 
black hole of mass equal to the mass coordinate. The thin long- 
dashed line gives the same but for a maximally rotating Kerr 
black hole. The thin dotted line denotes the specific angular 
momentum for the last stable orbit at the given mass of a black 
hole, assuming that all mass below forms a rotating black hole. 
Here, if the contained angular momentum is larger than that of 
a maximally rotating black hole, the black hole is assumed to 
rotate maximally. 



the core and the hydrogen envelope tends to render stronger 
differential rotation, coupling of the core with the envelope by 
magnetic torques enforces efficient transport of angular mo- 
mentum from the core. In Seq. A4, on the other hand, a de- 
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Fig. 6. Elapsed time from core helium exhaustion to the core 
carbon exhaustion, for Seqs. Al (20 M ), A2 (30 M Q ), A4 
(40 M ), A5 (50 M ), and A6 (60 M ). 

crease of the core angular momentum occurs mainly due to the 
stellar wind mass loss. In this sequence, significant differen- 
tial rotation only appears during the carbon core contraction. 
However, carbon burning starts soon after the central helium 
exhaustion (about 4000 yr after), and a large amount of angular 
momentum is still preserved in the core by the central carbon 
exhaustion. 

A comparison of the core angular momentum with the an- 
gular momentum of the last stable orbit around a black hole in 
Fig-Elreveals that the Seq. A4 fulfills the angular momentum 
constraints of the collapsar model. Table 1 indicates that indeed 
Sequences A2, A4, A5 and A6, as well as CI and C2 conserve 
enough core angular momentum to produce a collapsar. 

On the other hand, the core of the 20 M model of Seq. A2 
loses more angular momentum than the comparable more mas- 
sive models (see Table 0, and does not retain enough angular 
momentum to make a collapsar. This is mainly due to the rather 
long time scale of CO core contraction (Fig. |6j and the com- 
parably large envelope-to-core mass ratio in this case. This im- 
plies the possible existence of a lower mass limit for collapsar 
formation through the considered evolutionary channel, which 
may restrict collapsars to CO cores more massive than 10 M . 

The CO core mass of model Seq. A6 is about 40 M , which 
roughly coincides with the lower mass limit of pair instabil- 
ity supernovae (e.g. El Eid & Langer 1986 Kippenhahn & 
Weigert 1990 1. Except in a narrow mass range close to the 
mass limit, where pulsational pair instability may induce a 
large amount of mass loss, no stable iron cores are formed in 
the pair instability regime. Black holes may still form for very 
high masses (cf. Heger et al. 2003). However, for a normal ini- 
tial mass function, CO cores of about 40 M may be considered 
as upper mass limits for collapsar production via chemically 
homogeneous evolution. 

In Seq. B2, where metallicity is hundred times higher (Z = 
0.001) than in Seq. A4, quasi-chemically homogeneous evolu- 
tion still occurs, but the mass loss in this 40 M sequence spins 
down the core to angular momentum values below the collapsar 
requirement. Sequences C2 and C3, however, show that with 
the currently best estimate for helium star winds (Vink & de 
Koter 2005 1, collapsars are still being produced at Z = 0.001 
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Fig. 7. Specific angular momentum as a function of the mass 
coordinate, in the last models of Seqs. CI (upper panel) and C2 
(lower panel) [ after core carbon exhaustion in Seq. C 1 and near 
the core carbon exhaustion (Xq = 10~ 4 ) in Seq. C2]. The thick 
solid line denotes the mean specific angular momentum over 
the shells, as function of the mass coordinate. The thick dashed 
line gives the equatorial specific angular momentum. The thin 
two-dotted-dashed line denotes the angular momentum of the 
last stable orbit around a non-rotating Schwardzschild black 
hole of mass equal to the mass coordinate. The thin long- 
dashed line gives the same but for a maximally rotating Kerr 
black hole. The thin dotted line denotes the specific angular 
momentum of the last stable orbit at the given mass of a black 
hole, assuming that all mass below forms a rotating black hole. 
Here, if the contained angular momentum is larger than that of 
a maximally rotating black hole, the black hole is assumed to 
rotate maximally. 

(Table Fig. [7}. Given the still significant remaining uncer- 
tainties in Wolf-Rayet wind models, which partly imply very 
large clumping factors and thus further reductions in the Wolf- 
Rayet mass loss rates (Graefener & Hamann 2005 1, it appears 
conceivable that the considered evolutionary channel might be 
relevant for GRBs not only at high redshifts, but also in the 
local universe. 



4. Discussion 

Here, we investigate the evolution of rotating massive sin- 
gle stars at low metallicity, with the Tayler-Spruit dynamo 
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included. Our results show that if the initial spin rate is 
high enough, such that the time scale for rotationally in- 
duced mixing becomes shorter than the nuclear time (Maeder 
& Meynet 2000), the massive star may evolve in a quasi- 
chemically homogeneous way to become a rapidly rotating WR 
star on the main sequence (Maeder 1987 Langer fT992t . In par- 
ticular, it is shown that for low enough metallicity, this type of 
evolution can lead to retention of sufficient angular momentum 
in CO cores in the range 1OM . . .4OM to produce gamma- 
ray bursts according to the collapsar scenario (Woosley 1993 1. 

The critical values of metallicity and initial spin beyond 
which collapsars can be produced depend on uncertainties in 
the description of rotationally induced mixing and stellar wind 
mass loss. The Wolf-Rayet wind mass loss rate may be of cru- 
cial importance, as illustrated above. In general, more collap- 
sars are produced through chemically homogeneous evolution 
for lower metallicity, since less mass loss by stellar winds will 
keep the star more rapidly rotating, which favors mixing. 

However, due to the uncertainties in the stellar wind mass 
loss rates, and the uncertain distribution of initial stellar ro- 
tation rates (Maeder & Eenem2004i, chemically homoge- 
neous evolution at near solar metallicity can not be excluded. 
Observations actually provide evidence for homogeneous evo- 
lution in some young star clusters in the Galaxy and the 
Magellanic Clouds (e.g. Howarth & Smith 120011 Bouret et 
al. 120031 Walborn et al. 2004). Although such findings are in- 
frequent, this may in fact be consistent with the observation- 
ally implied low frequency of GRBs. Importantly, the results 
of Seq. CI and C2 with Z = 0.001 indicate that the blue stars 
as found by Bouret et al. (120031 and Walborn et al. (|5004 1 in the 
Magellanic Clouds which have the signature of homogeneous 
evolution might be promising progenitors of GRBs. Our solar 
metallicity Sequence B2 shows, on the other hand, the quasi- 
chemically homogeneous evolution does not always lead to a 
collapsar, in particular at high metallicity. 

It appears likely that at sufficiently low metallicity, the 
chemically homogeneous evolution channel will dominate over 
binary collapsar production channels as the helium star merger 
scenario or the tidal spin-up of helium stars — which are 
considered as a major path towards collapsars (e.g. Izzard et 
al. 1277041 Podsiadlowski et al. I2U041 Fryer & Heger 120051 : 
very metal poor stars may not evolve into red super-giants, and 
common envelope evolution without merging of the two stars, 
which is required in both mentioned scenarios, may be difficult. 
At zero metallicity, it is even questioned whether binary stars 
can form at all (Bromm & Larson 2004). Therefore, if GRBs 
occur from the first generations of stars, chemically homoge- 
neous evolution may provide their major formation path. 
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